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ABSTRACT 

We report the first results of a multi-epoch search for wide (separations greater 
than a few tens of AU), low- mass tertiary companions of a volume-limited sample 
of 118 known spectroscopic binaries within 30 pc of the Sun, using the 2MASS 
Point Source Catalog and follow-up observations with the KPNO and CTIO 4m 
telescopes. Note that this sample is not volume-complete but volume-limited, 
and, thus, there is incompleteness in our reported companion rates. We are 
sensitive to common proper motion companions with separations from roughly 
200 AU to 10,000 AU (~ 10" -> ~ 10'). From 77 sources followed-up to 
date, we recover 11 previously known tertiaries, three previously known candidate 
tertiaries, of which two are spectroscopically confirmed and one rejected, and 
three new candidates, of which two are confirmed and one rejected. This yields an 
estimated wide tertiary fraction of 19. 5l| 2 %. This observed fraction is consistent 
with predictions set out in star formation simulations where the fraction of wide, 
low-mass companions to spectroscopic binaries is >10%, and is roughly twice the 
wide companion rate of single stars. 
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Introduction 



Formation simulations have had a difficult time modeling the very close separations of 
many spectroscopic binaries. A m echanism is needed t o draw angular momentum away from 

Recent star formation simulations 



an already close pair of objects (IKiseleva et al. 



( ISterzik fc Durisenl 120031 ; iDelgado-Donate et al. 



1998) 



2004 



Umbreit et al.l 120051 ) show that one 



potential mechanism for the transfer of angular momentum is through three-body interac- 
tions. The third bodies used, in these cases, are cool dwarfs. Cool dwarfs are stellar and 
sub-stellar objects with spectral types > M and masses less than a few tenths of a solar mass. 
Interactions between loosely bound or totally unbound low-mass objects can dramatically 
tighten already-close orbits. These results predict that spectroscopic binary systems should 
have a larger fraction of wide cool companions than so-called 'single' stars. 



The dynamic interaction simulations of lSterzik &: Durisenl ( 20031) andlDelgado-Donate et al 



( 120041 ) produce some testable predictions. The simulations of IDelgado-Donate et al.l (120041 ) 
predict that, if a cool dwarf is found to be in a stable, >10 AU orbit, its primary is frequently 
(~ 75%) a tigh t spectroscopic binary. A similar qualitat ive result i s foun d in the dynamic 
simulations of ISterzik fc Durisenl ( 120031 ). Recent work by lLaw et al.l ( 120101 ) . which studied a 



sample of known, very wide M dwarf binary systems, found that 45l 16 % were higher order 
multiples. This supports the simulation predictions. Note that 'wide' in this work means 
separations > 10's of AU. 

The Delgado-Donate et al. simulations found that ~10% of their tight multiple sys- 
tems survive with wide, cool dwarf companions by the end of their simulations, 10.5 Myr. 
There are many e mpirical measurements of the cool dwarf companion frequency to stars. 
Gizis et al.l (120011 ) estimated a wide cool dw arf companion frequ ency of 18%±14%, based 
on only three L and T dwarf companions. ICarson et al.l ( 120091 ) examined 21 FGK stars 
within 20 pc, and found no co mpanions in a range of 20 AU - 250 AU down to masses of 



50 Mj. lLowrance et al.l (120051 ) used the Hubble Space Telescope to examine 45 young stars 



(~0.15 Gyr) at separations from 15 AU to 200 AU and masses dow n to typically 30 Mj to 
find or confirm 8 cool companions. McCarthy fc Zuckermanl ( 120041 ) examined a sample of 
~300 single, G, K, and M stars at separations from 75 AU to 1200 AU and masses as low 
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as ~5 Mj and found a cool dwarf binary frequency of 1-2%. iMetchev fc Hillenbrand! ( 120041 
20061 . 120091 ) examined ~250 'Solar Analogs' (F5-K5 stars). They found two new brown dwarf 



and 24 new stellar companions with separations from 28 AU to 1590 AU and probed masses 
down to ~10 Mj. They calc ulate an ultracool companion f requency of ~ 3%, which is 
marginally higher than that of McCarthy fc Zuckermanl ( 120041 ) . 



All of the above works are consistent with a low, wide, cool dwarf companion fraction 
of only a few percent. However, the primaries are al l apparently sing l e stars , which does not 
test the simulation results on tertiary companions. iTokovinin et al.l (120061 ) conducted such 
a study of 165 spectroscopic binaries, in which they searched for wide companions using the 
2MASS database. A subset of those objects (62) were observed at high spatial resolution 
with NACO on the VLT. They found a very high tertiary rate, adjusted for incompleteness, 
of 63%. They also found that the fraction of spectroscopic binaries with tertiary components 
is a strong function of spectroscopic binary period. Those with very short periods (less than 
12 days) almost all have wide companions (96%). 

It should be noted that other groups have conducted common proper motion (CPM) 
comparisons between various wide-field surveys, though none has specifically targeted wide 



companion s to spectroscopic binaries. The Brown Dwarf Kinematics Project ( IFaherty et al. 



20091 . |2010( ) has studied the kinematics of ultra cool dwarfs and found or co nfirmed several 

wide c ool companions to stellar primaries, including one spectroscopic binary. iLepine fc Bongiorno 
( 120071 ) performed a re-analysis of the Digitized Sky Survey using the custom built software 
package SUPERBLINK. They found that ~9.5% of Hipparcos stars have companions at sepa- 
rations wider than 1000 AU and proper motions great than O'.'15/yr. There have been several 
pa pers cross-corre l ating 2MAS S and SPSS which have reported individual discoveries, such 



MS 



Metchev et al.l (120081 ) and iGeiBler et al. 



( 201lh . T here have also been l a rge sy stematic 
cross-correl ations, such aslKirkpatrick et al.l ( 120101 ) and lSheppard fc Cushingl (120091 ). 2MASS 
to SDSS; or Deacon et al. (12009 ) . 2MASS to UKIDSS. However, these works focused on find- 
ing individual field objects wit h large proper mot ions, not finding multiple systems. Finally 
there is the Slowpokes survey (IDhital et al.ll2010l ) which cross-correlates USNO-B to SDSS 
and was designed to look for CPM companions. They focused on the general field popula- 
tion and only searched separations < 180" . This is narrower than our search radius and is 
limited to the optical and will not have the same sensitivity to extremely cool objects as our 
near-infrared survey. 

Here we report the first results from a study to measure the wide tertiary fraction 
around spectroscopic binaries via CPM, using 2 MASS as a first epoch and our own deep 
near-infrared imaging as the second epoch, which ITokovinin et al.l (120061 ) did not carry out. 
Section 2 describes the experimental setup and sample selection, while Section 3 outlines the 
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wide-field NIR imaging campaign and the data reduction procedures. Sections 4 and 5 detail 
our CPM analysis techniques and discusses the results, respectively. Section 6 summarizes 
the current work and our results. 



2. Experimental Design and Sample Selection 



We required our sample to differentiate statistically betwee n the 1 — 2% observed wide 



compani on rate of 'single' stars ([McCarthy Zuckerman 



dieted by lSterzik &: Durisenl (120031 ) , iDelgado-Donate et al 



2004) and the > 10% resu l ts pre - 



( 120041 ) . and lUmbreit et all ( 120051 ). 



Assuming Poisson statistics, a 1% — 2% frequency can be distinguished from a 10% frequency 
at the 3a level when ~100 objects are observed. Our spectroscopic binary sample of ~100 
was selected to detect bright to very faint wide companions using a CPM method. 

We used the 2MASS Point Source Catalog, which has astrometric uncertainties of ~0'/2 
and which was taken between 1997 and 2001, to provide the first ep och. Second epoch near- 
infrared imagin g was conducted w i th the KPNO 4m Flami ngos ( Elston et al.l 120031 ) and 
CTIO 4m ISPI (IProbst et al.ll2003l ; Ivan der Bliek et al.l 120041 ) instruments. These cameras 
have plate scales of O'.'3156/pix and O'.'3/pix respectively. Given our previous experience at 
measuring astrometry, we expected to be able to measure the position of objects in our 
fields to within ~10% of the typical seeing measurements. So, with the average seeing 
at KPNO and CTIO of ~ l'/O, we expected to obtain astrometry good to ~ O'/l. For a 
robust CPM measurement, we chose an approximately 5a shift in position to reject chance 
alignments with background objects whose motion closely matches that of the primary's. 
This is particularly relevant for those spectroscopic binaries whose motion is relatively small 
because our measurement uncertainties and systematics can make another, slower-moving 
background object appear to be a genuine CPM companion. A ~ 5a detection of CPM 
therefore requires a minimum motion of O'/l/yr. 

We also wanted our data to be sensitive to objects fainter than 2MASS by two to three 
magnitudes (similar sensitivities as UKIDSS lLawrence et al.l (120071 ) ). In this way we would 
be able to use the data reported here as a first epoch of deep observations sensitive to brown 
dwarfs considerably cooler than known T dwarfs and, pote ntially, the long sought after, and 



recen tly discovered, new spectral type range of Y dwarfs ( iKirkpatrickl l2008t ICushing et al. 



201 lh . As stated on the Flamingos Website, the ISPI exposure time calculator may be used 
for both instruments. The calculator indicates that a 15 minute exposure time yields 5a 
detection limits of 20.1, 19.5, and 18.8 at J, H, and K s respectively. These detection limits 
also represent significantly lower mass limits, for typical Galactic disk ages of a few Gyr, 
from our observations (~10 — 20 Mj) than is possible with 2MASS alone (~50 — 60 Mj). 



-5 - 



Also note that iFaherty et al.l ( 120091 ) has found that the typical proper motion of >M7 dwarfs 
with distances around 15 pc — 30 pc is >0'!2 /yr. This motion is well above the limit of our 
survey and the distance range corresponds to objects we could easily detect in our images. 
Thus, we will also be able to select for fast-moving objects in our fields to search for nearby 
very cool objects. 



The sample was selected from 'The Ninth Catalogue of Spectroscopic Binaries' (IPourbaix et al 



20041 ) (SB9). This catalog is a compilation of known binaries from the literature and comes 
with a few significant biases. When the SB9 catalog was first compiled (2004), there were 
about 1200 potential new systems in the literature that had not yet been evaluated for in- 
clusion. Additionally, the surveys for spectroscopic binaries th at form the core of th e SB9 



catalog are biased against fainter, lower-mass binary systems. IPourbaix et al.l ( l2004j ) indi- 
cates that this bias begins around K spectral types. Thus, while we treat this sample as a 
whole, it does not necessarily refle ct the statistics in a vo lume-limited sample which would 



be dominated by K and M dwarfs (IBochanski et al.ll201ll ). 



We ran the full catalogue of 2386 binaries through the Vizier interface of the Hipparcos 
catalog to determine the distance and proper motion of each object. Given the criteria 
derived at the beginning of this section, a volume-limit of 30 pc yielded an optimal sample 
size of 118 targets with appropriate motions (see Figured], left-hand panel), which are listed 
in Table [TJ By choosing nearby targets, we are sensitive to some of the brightest, lowest-mass 
brown dwarfs in the sky. Our data are able to find companions with separations greater than 
~10". This corresponds to projected separations greater than ~200 AU and periods greater 
than 1000 years for the typical primary distances and masses within our distance limit of 
30 pc. 

The resultant sample is both volume and proper-motion limited to those sources in the 
Hipparcos database. This introduces some biases, including kinematics: only objects moving 
faster than ~ 14 km/s at 30 pc were selected. This kinematic bias favors the selection of older, 
faster moving systems. Thus, we are biased against detecting the youngest systems with 
the lowest-mass companions, as they will have cooled to temperatures below our detection 
limit. However, without knowing a priori what the ratio of young systems to old systems 
is, this bias is very difficult to estimate. There is also a spectral type bias against M dwarf 
primaries, due to the V magnitude limit of Hipparcos and the relative lack of data on M 
dwarf spectroscopic binaries. Also, note that this sample is not volume-complete. As can 
be seen in the right hand panel of Figure [JJ the number of objects continues to increase to 
distances of roughly 20 pc and then levels off. This is likely a result of the magnitude limit of 
the Hipparcos catalog. We are biased again the fainter, lower mass, more distant primaries. 
These biases will be quantitatively accounted for via a Bayesian/Monte Carlo simulation in 
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a future paper. 

Figure [2] displays two histograms of the spectral type distribution of our primary sample. 
The left-hand panel of Figure [2] shows the overall spectral type distribution, and we can see 
that we have mostly F, G, and K stars. Thus, our sample will most accurately examine 
the tertiary fraction of spectroscopic binaries composed of F, G, and K stars. The right- 
hand panel of Figure [2J shows how the spectral type distribution chang function of 
distance, which demonstrates that our sample is magnitude limited. This makes sense, given 
the selection biases of the SB9 catalog towards brighter, more distant, earlier spectral types. 
Thus, within our sample, we are likely missing G and K-type spectroscopic binaries in our 
larger distance bins. Given that we have a comparable overall number of F, G, and K stars 
in the sample, we are probing the statistical companion rate as a function of spectral type 
equally across those types. 

Another factor to be taken into consideration is that of orbital motion. We expect 
to locate companions as close as 10" — 20" from the central binaries. Given our smallest 
target distance of around 5 pc, the minimum angular separation of 10", and the typical time 
baselines of 10 years, the orbital motion is maximal at ~120 mas/yr. This is on the same 
order as our typical 2o CPM search criterion limit we derive in Section H~T1 of ~ 100 mas/yr. 
However, we only have a handful of systems at distances less than ~10 pc that would have 
such large maximal orbital motions. This will be discussed in more detail in Section I4.1[ but 
we conclude that the orbital motion effects will be minimal. 



3. Observations and Data Reduction 
3.1. Observations 

The observations described here were carried out during four successful observing runs: 
May 3-6 2007 (KPNO); June 3-6 2007 (CTIO); January 20-23 2008 (CTIO); January 28-30 
2008 (KPNO). Table [1] lists when and where each target was observed. 

We obtained standard calibration data in the afternoon before each night's observing, 
including sets of dome flats in all three bands (JHK S ), bias frames, and dark frames. Our 
science data consisted of dithered JHK S imaging of each target field. Dithering allowed 
the construction of accurate maps of the IR background in each exposure which were then 
subtracted from each individual frame. We set individual dither position exposure times to 
keep the background well below the beginning of the non-linear regime (~10000-13000 ADU 
for both instruments). We could not avoid saturating on the very bright central spectroscopic 
binaries. The saturated areas typically effect the inner part of the image to radii of 10 - 15 
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arcseconds but could reach to much larger radii (see section l5TTj) . 

Using IPSI, we modified a 15-point random dither script for each band pass and adjusted 
the exposure time and number of coadds according to the background at that band pass each 
night, such that the total exposure time per pointing was 60 seconds. For Flamingos the 
procedure was slightly different, due to the lack of a coadd feature. We instead repeated our 
dither patterns, 5X5 for the May 2007 run and 4X4 for the January 2008 run. We selected 
appropriate exposure times to keep the background in the linear regime. This required us 
to repeat the dither pattern two to four times. 



3.2. Data Reduction 

3.2.1. Basic Reduction 



Data reduction procedures for both instruments were performed identically, using IRAF 
tools developed for ISPlE and the publically available software packages WCSTOOLS ( Mink 



20021 ) and SWarp. For each night of data, we first created a bad pixel mask, usually from a 
dome flat image, using the IRAF task ccdmask. This mask was applied to all data for that 
night, including calibration data, using the IRAF task fixpix. This task linearly interpolates 
over the marked bad pixels from ccdmask using nearby 'good' pixels. We reduced the JHK S 
data for each science target separately in the following manner. First, we dark-subtracted 
the raw frames and corrected for non-linearity. Next, using the xdimsum IRAF package, we 
estimated the sky background, subtracted it, and masked the holes (the result of subtracting 
dithered images with stars in them from one another). In the final step, we flat-fielded the 
sky and dark-subtracted images. 



3.2.2. Mosaiking and Coordinate Solutions 

Since the main goal of this project was to obtain accurate proper motions for the objects 
detected in our fields, we needed to combine the individual frames into a final image and refine 
the resultant coordinate solution carefully. We first oriented each final reduced image, such 
that North is up and East is left. Sources were then selected to create a rough coordinate 
solution. This was done using the starfind IRAF task. We selected these initial lists to 
be unsaturated sources in each field of medium brightness, typically from 10 — > 14 mag 



3 http:/ /www. ctio.noao.edu/instruments/irjnstruments/ispi/New/UsersGuide/datared. html 
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at J-band, to match against the 2MASS PSC from the full release (ISkrutskie et al.ll2006l ). 
We then fed the output from starfind into the WCSTOOLS task imwcs, which provided a 
preliminary coordinate solution. 

Next, we refined the rough coordinate solution provided by imwcs, and corrected dis- 
tortion across the field of view using the IRAF task ccmap, which we ran interactively for all 
images via the cirred package task do_ccmap. This was done to determine the average dis- 
tortion solutions across all of the images. In ccmap, distortion was quantified by polynomial 
fit in x and y pixel positions, with orders set to either four, for relatively sparse fields (those 
with fewer than ~40 sources, which corresponds to roughly a third of our observed sample), 
and six, for more populated fields. When ccmap was run on each image, we interactively 
removed sources that have large differences from their 2MASS coordinates until the residuals 
of the coordinate solutions were on the order of O'.'l. The number of sources removed varied 
from image to image and could be as few as 4 or 5 or as many as 30 or 40 for very crowded 
fields. This was never more that ~10% of the sources in any given field. The removed 
sources were typically mismatched artifacts due to the regions around the bright primary. 
Note that the residuals we obtained were similar to those of the 2MASS PSC and matched 
our expected astrometric precision. Also, similar coordinate solutions and final coordinates 
were found in all three bands, J, H, and K s . Finally, we used SWarp, including the average 
distortion correction created via do-ccmap, to perform the final mosaiking. 



4. Identification of CPM Companions 
4.1. Technique 

We measured our second epoch astrometry from the J-band mosaicked images for the 
majority of our sources. Again, we selected sources using the IRAF task starfind, extracted 
photometry using phot, and measured RA and DEC coordinates from the image WCS in the 
fits header, using wcsctran. We manually cleaned this catalog of sources of all spurious ob- 
jects by visual inspection, using ds9. The majority of these spurious sources were associated 
with the bright halo of the saturated central spectroscopic binary and typically numbered 
between a few dozen and several hundred per field. We then compa red the cleaned list o f 



sources to objects found in the field selected from the 2MASS PSC (ISkrutskie et al.ll2006l ) 
using custom-built Fortran scripts. 

Comparison of our astrometry to that from 2MASS was a two-fold process. The first run 
looked for objects that did not move (defined as objects in our data that are within 0'.'5 of their 
2MASS positions) and marked them as matched. This helped to prevent mismatches with 
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nearby faint sources when looking for moving objects. The second run ignored these matched, 
non-moving objects and looked for objects with motions of a similar magnitude to that of 
the target spectroscopic binary in the field. Our program then filtered the list of 'moving' 
objects further by determining the uncertainties on all the proper motion measurements for 
that field. The majority of objects in a given field display little to no motion, and those 
measurements follow a normal distribution. Thus, we can use the standard deviation of that 
normal distribution as an estimate of the uncertainty of our proper motion measurement. We 
used the 2a limit as the threshold to determine the proper motion uncertainty in both \x a and 
Us in mas/yr. The limits ran from 62 mas/yr up to 160 mas/yr with an average value of ~100 
mas/yr. These limits naturally include the ~ (X'l residuals from the coordinate solution. We 
considered any object whose motion falls in that 2a to be a CPM companion candidate. 
Example plots of the CPM diagrams for fields both with and without a companion can be 
seen in Figure El 

As mentioned in Section 2, the orbital motion of candidate companions must be exam- 
ined. We determined the maximum orbital motion of a companion in all of our fields by 
taking the typical inner working angle of ~10", the distance to each primary, and typical 
primary masses of ~2 M and applying Kepler's 3rd Law. These maximal values were then 
compared to the 2a proper motion uncertainty values. In all cases, save one, the maximum 
orbital motion was less than the uncertainty used to select candidates. Only in the field 
of HIP 88601, our closest observed target, was the maximal motion of 123 mas/yr greater 
the proper motion uncertainty of 120 mas/yr. Additionally, several of our targets that fall 
within 10 pc also have maximal orbital motions that are within a factor of 2 of the proper 
motion uncertainties. However, at radii of 15" to 20" the orbital motion falls well below the 
uncertainties. These fields represent about 10% of our observed sample, but it is only a rel- 
atively small fraction of each field in which we could have thrown out a potential candidate. 
Thus, we most likely did not remove any candidates because of their orbital motion. 

Despite the above selection thresholds, mismatches can still occur, which are typified 
in the left-hand panel of Figure [3j To remove this last set of spurious sources, we employ 
one further candidate selection criterion, visual comparison of the Flamingos/ISPI images 
with the Digitized POSS plates. Given the very red optical-NIR colors of brown dwarfs, we 
do not expect them to have counterparts in the POSS R-band. Nearly all of the candidates 
that passed the earlier criteria we found to be spurious. They typically occurred in fields 
around primaries with higher proper motions (> ~0"4/yr), as a result of mismatches between 
brighter 2MASS sources with small motions (hence not thrown out in the first sweep) and 
faint sources from our Flamingos/ISPI images. These non- moving, faint sources in our data 
were matched to brighter 2MASS sources because they are too faint to appear in the 2MASS 
PSC. This leads the matching software to find the closest available unmatched 2MASS source, 
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which is always a nearby bright source. These sources are also typically matched against 
a faint source in the Digitized POSS, and the bright 2MASS source has a corresponding 
bright source in our Flamingos/ISPI data. So, we did not discard any potentially very red 
companions. 

Through the visual comparison of each candidate to the DPOSS R-band plates we 
reduced the list of candidates from dozens to 13, which are listed in Table [2] along with 
five other known bright companions that are saturated in our data and thus not detected 
in this analysis. We uncovered the five saturated companions via a literature search using 
the SIMBAD database for each of the target spectroscopic binaries we observed. These 
companions are all very bright, with apparent J magnitudes less than 8. For those known 
objects whose motion we measured, we successfully recovered their known proper motions 
within our uncertainties. This was a final verification step for the coordinate solution process 
discussed in Section 13.2.21 



4.2. Candidate Follow-up 



Among our CPM companions, we identified 14 that had been previously noted as ver- 
ified or candidate companions in the literature. Of these, 11 are confirmed companions, 
discovery references are given in Table [2J based on their motions and spectral types, the lat- 
ter providing a spectrophotometric distance that is consistent with the Hipparcos distance 
of the spectro scopic binary. Three additional ca ndidate companions were also recovered 
(HIP 41211C JTokovinin et al.l boOfih: HIP 83608C jPommanget fc NvsH2002h : HIP 101769C 
(iDommanget fc Nvslbool vh for which CPM or spectroscopic confirmation had not yet been 
obtained. Finally, we identified three new candidate companions. 

Candidates were followed up w ith low-resolution 1—2.5 /im spectroscopy using the prism 
setting of Spex feavner et aDl2003h on NASA's InfraRed Telescope Facility (IRTF). In all of 
the observations the spectrograph slit was aligned with the parallactic angle to minimize the 
effects of differential refraction. The data were taken in the ABBA dither pattern format to 
allow for easy subtra ction of the infrared background. We reduced these data using Spextool 
(ICushing et al.l 120041 ). an IDL based software package designed specifically for SpeX data. 
This reduction package removes the background, traces and extracts the spectral data, and 
wavelength and flux-calibrates the data with an AOV standard star. The final reduced and 
calibrated spectra are then compared to spectral standards from the IRTF Spectral Library 
(ICushing et al.ll2005t iRayner et al.ll2009l ). and are displayed in Figures 0] - EE 
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4-2.1. Previously Identified Candidate Companions 



HIP 4 121 1C: The closer com panion to HIP 41211, in Tabled was noted as a visual compan- 
ion in lTokovinin et al.l ( 120061 ). Our data demonstrated CPM, and our SpeX prism spectrum 



yielded a spectral type of M4.5 (Figure H]), which is consistent with being a physical compan- 
ion. This c o nfirm ation, in addition to the very wide M5.5 companion to HIP 41211 found by 
Reid et al.l ( 120031 ). makes this system the only confirmed quadruple currently known in this 
sample. The projected separation for the newly confirmed companion is 980 AU, whereas 
the companion found in iReid et al.l (120031 ) is nearly ten times as large, 9640 AU. This sys- 
tem appears to be hierarc hical and is a stable configuration according to the criteria of 
Eggleton fc Kiseleval (119951 ) . The primary spectroscopic binary in this system has a metal- 
licity measu rement of \Fe/H] = — 0.31 relative to the Sun, along with an age estimate of 



3.41^ Gyr (iNordstrom et al.ll2004f ). 



HIP 83608C: The companion to HIP 83608 was noted as a visual c ompanion in the Catalogu e 
of the Components of Double and Multiple Systems (CCDM2) (jPommanget &: Nysl 120021 ). 
Using our SpeX spectrum of this object, we found that the spectral type is definitely that 
of a M dwarf but that it is not a good fit to a single spectral type (see Figure [5]). The bluest 
part of the spectrum is better fit with a M4 dwarf, whereas the bulk of the spectrum mor e 
closely resembles a Ml dwarf. Using the Mj-SpT relationship from lHawley et al.l (120021 ). 
we estimated the expected spectral type, given an absolute J magnitude of 6.9±0.04, as 
M1.5-M2. This agreed well with the overall Ml spectral type estimate. However, there was 
still the discrepancy in the blue portion of the spectrum. 

We eliminated the possibility of this object being a background, low-metallicity subdwarf 
by placing it on a red uced proper motion diagram. We used the prescription laid out in 
Lepine fc Sharal (120051 ) for Hy versus V — J where 



H v = V + 5 * logM"/yr)) + 5. 

HIP 83608C has a V magnitude in the CCDM2 of 13.8; so, with a V-J color of 4.78 and 
a total proper motion of / . / 1164 / yr±0 / / 043/yr, we got a Hy of 14.1. When this was placed 
on Figure 30 of iLepine fc Sharal (120051 ). it lies just above the galactic disk dwarf sequence 
by about a magnitude. Thus, it is not a chance background subdwarf, and, given its slight 
over-brightness, it could be a binary itself. If a binary, this system would likely be composed 
of an early M dwarf, given the majority of its spectrum, and a later type M dwarf as the 
projected absolute J magnitude matches well with the estimated spectral type. Note that 
our Flamingos data do not show any elongation of the psf. More data, particularly an 
optical spectrum, are required to determine the full nature of this companion. The primary 
spectroscopic binary in this system has a metallicity of [Fe/H] = —0.01 relative to the Sun 
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along with an age estimate of 2.2±g;f Gyr (INordstrom et all 12004 ) . 

HIP 101769C: The CPM object in the field of HIP 101769 was noted as a visual companion 
in CCDM2 and is known as CCDM J20375+1436C. The follow-up SpeX spectrum, displayed 
in Figure [6], shows that this object is a G dwarf and is more likely to be a chance background 
alignment than a physical companion. The probability of a given object being a background 
interloper is discussed further in Section [5j 



4-2.2. New Tertiary Candidate Companions 



HIP 72603C: The companion candidate to HIP 72603 has near-IR colors and absolute mag- 
nitudes consistent with an early-T dwarf. However, it is also clearly detected in the DPOSS 
BRI plates. Its SpeX prism spectrum most closely matches an early M dwarf standard (Fig- 
ure [7]). We also performed the same reduced proper motion analysis as for HIP 83608C. We 
found that this object has a H v = 16.8 and a V — J = 3.1. When placed on Figure 30 of 
Lepine &: Shard (120051 ) . it falls in the subdwarf regime. Thus, it is not a physical companion 
but a chance background alignment. 

HIP 86722C: The new tertiary companion candidate to HIP 86722 has near-IR absolute 
magnitudes and colors consistent with a mid-M dwarf at the same distance as its putative 
primary. It is also detected in DPOSS. The candidate displays consistent CPM measurements 
over 50 years, including data from both POSS-I (1950s) and POSS-II (early 1990s). We 
obtained a SpeX prism spectrum and determined its spectral type to be M4.5 (Figure EJ). 
Thus, this object appears to be a genuine physical companion. The primary spectroscopic 
binary in this system h as a metallicity measurement of [Fe/H] = —0.39 relative to the Sun 
(INordstrom et al.ll2004h . 



HIP 97944C : The new tertiary companion candidate to HIP 97944 has near-IR absolute 
magnitudes and colors consistent with a mid-M dwarf at the same distance as its putative 
primary. As with HIP 86722C, this candidate is also detected in DPOSS, in both POSS- 
I and POSS-II, and displays consistent CPM over a baseline of 50 year s. Optical I-band 
photometry {mi = 11.94 ± 0.03 mag) from DENIS (jEpchtein et al.lll997l ) and the resultant 
optical-NIR colors are also consistent with a mid-M dwarf. We obtained a SpeX prism 
spectrum and determined its spectral type to be M5 (Figure [HD • Thus, we find this object 
to be a genuine physical companion. The primary spectroscopic binary in this system h as a 
metallicity measurement of [Fe/H] = —0.11 relative to the Sun (jNordstrom et al~ 2004 ). 
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5. Discussion 
5.1. Sensitivity 

We made preliminary estimates of our survey sensitivity, both in separation from the 
bright spectroscopic binary primary and in magnitude. This analysis was performed by 
creating a psf star for each field. This psf star was then randomly placed at ten locations 
within the field, each with a random magnitude. The daophot IRAF package was used to 
create the psf star and place the random, 'fake' stars. The resultant images were searched 
for point sources using identical starfind parameters from the initial search for candidates, 
as described in Section 14.11 The insertion of ten fake sources at a time was repeated 1000 
times for each field for a total of 10,000 fake sources. 

The magnitude limit was then determined as a function of separation from the central 
binary by comparing the number of fake sources inserted to the number recovered by our 
search parameters. The left-hand panel of Figure [10] displays the sensitivity curve generated 
for the field of HIP 39064 at 50% and 90% completeness levels. For this field, we see that 
the inner 10-15 arcseconds is mostly lost, due to the bright primary. However, outside of 
~15 arcseconds, our sensitivity is fairly uniform, with an average 50% completeness of J 
= 18.3 mag and a 90% completeness of J = 17.6 mag. In comparison to the field of HIP 
39064, which corresponds to the typical brightness of our central binaries (V = 7.70 mag), 
the sensitivity curve of HIP 44248 (V = 3.97 mag) is displayed in the right hand panel of 
Figure [TU1 This is one of the brightest objects in our sample. The average 50% completeness 
is 18.1 mag and the average 90% completeness is 16.7 mag. However, as expected for a 
brighter central object, the radius at which a uniform sensitivity is reached is much larger, 
30 — 40 arcseconds. 

It should also be noted that the effects of the bright central binaries are not limited to 
the central 10-15 arcseconds. Among the brighter primaries, we noted a ringing effect in our 
images. These looked like ripples in a pond, centered on the primary. This ringing also caused 
wide swings in sensitivity across the FOV. This effect can be seen in the sensitivity curve of 
HIP 44248 in the right-hand panel of Figure [101 It is particularly noticeable between 20 and 
50 arcsecond separations where the 90% completeness limit jumps by several magnitudes a 
couple of times. While we are not completely certain what the cause of the 'ringing' is, it is 
likely due to scattered light within the camera from the extremely bright central objects in 
our fields. 

Figure [TT] displays a histogram of the 50% and 90% completeness limits for all fields 
observed in this work. The median 50% limit is 19.1 mag and the median 90% limit is 18.2 
mag. From these data, we can see that we did not achieve our expected sensitivity in many 
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fields (J band limit of 20 mag), but can consistently recover objects one to two magnitudes 
brighter. These sensitivity curves as a function of radial separation from the central binary 
will be used in a future statistical analysis of the sample. 



5.2. Overall Wide Companion Rate 



Our primary goal in this initial study was to test whether known spectroscopic binaries 
have an enhanced, wide tertiary rate, compared to that of other stars. Table |3] lists the 
frequency of wide tertiary companions in our sample broken down by spectral type of the 
primary member of the spectroscopic binary system. With the exception of the G and B 
stars, of which we only observed 4, the typical wide tertiary fraction is ~ 20 — 25%. The 
overall rate, which we determined by assuming that these systems are drawn from a binomial 
distribution, is 19.5± |j%. This is a prelimina ry rate, as we do not account for incompleteness, 
selection effects, etc. iTokovinin et al.l ( 120061 ) probes the same overall spectral type primaries 



as we do, but also to small er separations. In orde r to compare consistent samples we only 
select companions found in ITokovinin et al.l (120061 ) that we could have found in our sample. 
Since we probe the same primaries we only selected companions in the same separation range 
(> 10"), or 27 companions. This yields a tertiary fraction of 16.8 ± 2.9% from a sample of 
165 spectroscopic binaries, which is comparable to our findings. 

Our wide companion fraction can also be compared to that of apparently single stars for 
a similar range of primary spectral type. We also compared our com panion rate to that of ap- 
parently single stars. For example, iMcCarthy fc Zuckermanl (120041 ) found a wide companion 
rate of ~11%-12%, which is still lower than either our own results or those of iTokovinin et al. 



( 120061 ). Thus, both this work and that of ITokovinin et al.l ( 120061 ) are consistent in finding 
an enhanced tertiary companion fraction to spectroscopic bi naries with respe ct to 'single' 
stars. This trend continues into the M dwarf regime as well. lLaw et al.l ( 120101 ) surveyed 36 
known wide M dwarf binary systems to look for close companions to either members. They 
find that 45^}g% of their wide bina ries are actually high-ord er multiple systems. Thi s 



agree s 



wit h the simulation pred ictions of ISterzik fc Durisen d2003h. IPekado-Donate et all (120041) 
and lUmbreit et al.l ( 120051 ). as well as with our results. IMcCarthy fc Zuckermanl (j2004j ) sur- 
vey nearby solar type stars (FGK) over a wide range of separations and find a substellar 
companion rate of ~ 1%. This is similar, within the un c ertain ties, to our observed substellar 
companion rate of 3.8 ± 2.2%. iMetchev fc Hillenbrand! (120091 ) perform an exhaustive survey 
of 266 FGK stars with both AO imaging and wide field studies and find a substellar com- 
panion rate of 3. 2+2^% for separations of up to ~1600 AU. This rate is statistically identical 
to the rate we found. 
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The substellar companion fraction we report is most likely a lower limit, as we should 
find more faint companions when we obtain our second epoch deep imaging. Those data will 
have comparable sensitivities to our first epoch (J; im ~18.2 at 90% completeness for most 
fields outside of 10" to 20" separations) and will be sensitive to very faint T and later dwarfs. 
This is because we rely on 2MASS for ou r first epoch astrometry. 2MASS is fairly complete 
to distances of 25 - 30 pc for L dwarfs (ICruz et al.l 120071 ). However, T dwarfs, which all 
have absolute J magnitudes of ~ 14 — 16.5, are only detectable to distances of 15 - 20 pc, 
given the J-band 2MASS limit of ~16. So, at this time, it is difficult to tell if the substellar 
companion rate is enhanced. 

Finally we examined possible correlation between the mass of the primaries and the 
mass of the secondaries by using spectral type as a proxy for mass. There was no obvious 
correlation, although the majority of the tertiary companions were of the M spectral type 
for all primary spectral types. This was not true for G type primaries, whose only confirmed 
companions were L dwarfs. In any case, the large majority of wide tertiary companions 
noted in this wor k were of a considerably lower mass than their primaries, w hich fits with 
the predictions of ISterzik fc Durisenl (120031 ) and iDelgado-Donate et al.l (120041 ) . 



5.3. Background Interlopers 



We found two background interlopers in our CPM sample, one distant G dwarf (HIP 
101769C) and one background M subdwarf (HIP 72603C). HIP 101769 has a proper motion 
near the limit of our minimum motion of O'.'l/yr, while HIP 72603 has a proper motion much 
higher (~0'/2/yr). It has been found that the numbe r density of moving ob jects rises as the 
inverse cube of the magnitude of the proper motion (ILepine fc Shard 120051 ) . Thus, it is not 
surprising that the interlopers we found lie at the lower end of our motion spectrum. 



Further work by ILepine fc Bongiornol (120071 ) provides a quantitative mechanism for de- 
termining the likelihood that a given CPM candidate is a chance alignment of a background 
object when the overall motion of the object is > CfUb/yr. This analysis was based on the 
fact that the number of objects at a given proper motion increases with smaller motions 
and that the chance of a random alignment increases with greater angular separation. They 
derived the following formula to quantify these correlations: 

AX = [O/O.15)~ 3 - 8 A0A/^ 

where \i is the magnitude of the proper motion of the primary in arcseconds per year, A9 is 
the difference in angular position between the primary and the candidate CPM companion in 
arcseconds, and A/x is the difference in proper motion in arcseconds per year. The quantity 
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AX measures the likelihood that a given object is genuine. iLepine fc Bongiornol (120071 ) found 
that, when the value of AX is around 1, there is a 50% chance of the candidate companion 
being a chance alignment. This increases to well above 90% for values of AX > 1.2. 

The value of AX for the two interlopers we find in our sample are 1.2 for HIP 101769C 
a nd 3.9 for HIP 72603C. Sinc e the motion of HIP 101769C is below the limit of the analysis 



in 



Lepine fc Bongiornol (120071 ). it is not clear how effectively this formula can be applied. All 



of our other candidates have values under 1. 



5.4. Spectroscopic Binary Periods 



The comparison between the target sample of this work and that of iTokovinin et al. 



(120061 ) produced some interesting results, particularly when comparing the orbital periods of 
the target spectroscopic binaries. Tokovinin's program set out with the same goal as ours: 
examination of the tertiary fraction of spectroscopic binaries as a means of testing binary star 
formation simulations. They wanted to maximize the chances of detecting tertiaries; so, they 
selected only spectr oscopic binaries that only have periods of less than 30 days . The simu- 
lation predictions of ISterzik &: Durisenl (120031 ) and iDelgado-Donate et al.l (120041 ) argue that 
the Kozai mechanism can tighten these systems by transferring angular momentum from the 
tight system to a wide third member. Thus, the tighter systems should, pr eferentially, have 



a high er tertiary companion rate than wider spectroscopic binaries. In the ITokovinin et al. 



(120061 ) study they do find that tighter spectroscopic binaries tend to have a higher fraction of 
tertiaries and that the tertiary fraction rises to nearly 100% for spectroscopic binaries with 
periods of less than a day. 

Our selected sample is volume and minimum proper motion limited, while the period 
of the spectroscopic binary is unconstrained. Table H] lists companion fraction as a function 
of the period of the spectroscopic binaries in the observed sample, while Figure [12] displays 
a histogram of the data in Table HI It is broken down in log period bins centered on the 
values listed with ±0.5 dex widths. We calculated the fractions, assuming that the data are 
binomially distributed (Table H|). Our number of primari es as a function o f bina ry period 
is fairly flat from 1 day to 10,000 days. The sample from ITokovinin et al.l (120061 ) does not 
examine spectroscopic binaries with periods longer than 30 days, which corresponds to the 
last three bins of Table H] and Figure H21 The tertiary companion rate, however, peaks at 
the smallest period bin for which we have significant data, 50% ± 13.3% for spectroscopic 
binaries with periods between 1 and 10 days. The companion fraction then decreases for 
spectroscopic binaries with periods between 10 and 100 days, but then increases again to 
between 15% and 30%. Note that the most significant of these variations from our baseline 
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19.5% companion rate is the 50% r ate at small peri o ds, w hich is a ~2<r event. Thus, we 
can marginally confirm the result of iTokovinin et all (120061 ) . which found that the tertiary 
companion rate drops by about a factor of two from very close binaries (periods less than 7 
days) to wide binaries (periods between 7 and 30 days). A larger sample size is needed to 
provide a more robust measure ment of the varia t ion of the companion fraction as a function 
of spectroscopic binary period. ITokovinin et all ( 120061 ) also find that there is no correlation 
between the period of the binary and the separation of the tertiary, which we find as well 
(see Figure IT3"j) . However, our result of an increase in companion fraction for spectroscopic 



binari es with periods longer than 100 days is quite different from that of ITokovinin et al. 
(bood ). 



We found a significant fraction of wide spectroscopic binaries have tertiary companions. 
This demonstrates that enhanced wide companion rates also apply to wider spectroscopic 
binaries (Table HJ). However, this result is preliminary, as there is a large area of separation 
space that this work has not yet explored , particularly ter t iary c ompanions with separations 
less than 10". It should be noted that ITokovinin et al.l (120061 ) found nearly half of their 
tertiary companions in this separation range. 



6. Summary 



We surveyed a volume-limited sample of 77 spectroscopic binaries for CPM tertiary 
companions. We found or confirmed a total of 4 tertiary companions, two of which are 
brand new to this work. When combined with the 11 previously known companions, this 
yields a total of 15 companions in 13 triple systems and 1 quadruple (HIP 41211), for 
an observed wide companion rate of 19.5%+g^. Note that these numbers are valid for 
separations > 10 — 20 arcseconds or hundreds of AU. Also, our magnitude limit is up to 
3 — 4 magnitudes deeper than 2MASS; thus, we are currently sensitive only to sources also 
detected in 2 MASS. 

Initial comparisons with other observed tertiary or wide binary rates yie ld contradictory 



result s. First, when we compared these data to the only other similar survey, ITokovinin et al. 



( 120061 ). we found a similar tertiary companion fraction. When we compared these fractions 
to 'single' FGK stars, of which the bulk of our sample binaries are composed, we foun d an en- 
hanced fraction of wide companions (~11%-12% from McCarthy fc Zuckermanl (120041 ). versus 
19.5% from our work). However, whe n comparing the substellar ter t iary rate of this sample, 
3.8=fc 2.2%, to that of similar surveys ([McCarthy &: Zuckermanll2004l ; iMetchev &: Hillenbrand 
20091 ). 3.2%, we found little difference, though it should be noted that all cases deal with 
small-number statistics. These results are preliminary and will require additional analy- 
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sis to confirm. In particular, a second epoch of deep imaging will be ne eded to search for 



fainter tertiary com panions, which include objects in the Y dwarf regime (iKirkpatrick! 12008 



Cushing et al.ll2Qllh . 
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Fig. 1. — Left: Measured distances versus proper motions for all Hipparcos spectroscopic 
binaries within 30 pc of the Sun. Open triangles represent each object in the sample. The 
horizontal dotted line gives the minimum proper motion cut-off for inclusion in the final tar- 
get list. Right: Cumulative histogram of the number of spectroscopic binaries as a function 
of distance. Note that the number flattens at distances larger than 20 pc, indicating that 
this sample is volume-limited and not volume-complete. 
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Fig. 2. — Histograms of the spectroscopic b inary spectral types in the sample we derived 
from the SB9 catalog (IPourbaix et al.l 12004 ). The bottom panel shows our entire sample, 
note that we have a large number of FGK stars but very few B, A, or M stars, and no O 
stars. The top three panels show histograms of the relative fraction of primary spectral types 
in our sample broken down into three distance bins: pc - 10 pc (solid red); 10 pc - 20 pc 
(dashed green); 20 pc - 30 pc (dot-dashed blue). Note that the peak of the distribution shifts 
to earlier spectral types with increasing distance. 
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Fig. 3. — Representative plots of the CPM measurements for the fields around two target 
spectroscopic binaries, x's represent the measured motion of each object in the field detected 
in both the 2MASS PSC and in the ISPI/Flamingos data. The five pointed star indicates 
the proper motion of the primary, and the two large circles are the la and 2cr detection 
thresholds for candidate co-moving companions measured in mas/yr. The left panel is that 
of the field around HIP 75312 and it contains a known L dwarf tertiary. The right panel is 
that of HIP 105312 and is an example of a field without a companion. The HIP 105312 field 
also typifies the mismatches that lead to spurious high proper motion candidates discussed 
in Section 14.11 Note the source that is not clustered around zero motion; it is a spurious 
detection. However, since it is not within la — 2a of the primary it is not counted in this 
case. 
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Fig. 4. — SpeX prism spectrum of the candidate companion to HIP 41211 (black solid 
l ine) compared to M3-M6 dwarf spectra from the IRTF Spectral Library (red dotted line) 
(jCushing et al.ll2005l ; iRayner et al.ll2009l ). We find that the candidate is a M4.5 dwarf, based 
on this spectrum. Thus, coupled with the measured magnitudes and colors (see Table [2]), it 
is consistent with being a bona fide companion. 
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Fig. 5. — SpeX prism spectrum of the candidate companion to HIP 83608 (black solid line) 
compared to standard spectra (red dotted line), M0-M3 in the left pan el and M3-M6 in the 



right panel. We ob tained the standards from the IRTF Spectral library (ICushing et al.ll2005 
Rayner et al.l 120091 ). The overall morphology of the spectrum is that of a Ml dwarf in the 



red portion of the spectrum and a M4 dwarf in the blue. It could be a composite of an early 
and a mid M dwarf and, thus, a binary itself. 
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Fig. 6. — SpeX prism spectrum of the candidate companion to HIP 101 769 (black solid line) 



compared to stand ard stellar spectra from the IRTF Spectral Library ( jCushing et al.l 12005 



Rayner et al.ll2009l ) (red dotted lines). This object is clearly a late-G dwarf and, thus, cannot 



be a bona fide companion, given its projected absolute magnitudes and colors (Table [2]). 
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Fig. 7. — SpeX prism spectrum of the candidate companion to HIP 72603 (black solid 
line), com pared to M0-M3 dwarf spectra (red dott ed lines) from the IRTF Spectral Library 
standards (jCushing et al.ll2005l ; iRayner et al.ll2009l ). The overall morphology of the spectrum 
is that of a M dwarf, but it is not a good fit to any particular template. Coupled with 
the reduced proper motion determination (Section I5.3I) . it is likely a unrelated background 
sub dwarf. 



-29 - 



c 4 
o 

CO 

c 

o 
o 

+ 



a. 

CN 



X 

zs 



X 

zs 








1 .5 2 
Wavelength(/xm) 



2.5 



Fig. 8. — SpeX prism spectrum of the candidate companion to HIP 86722 (black solid 
l ine) compared to M3-M6 dwarf spectra (red dotted lines) from the IRTF Spectral Library 
(ICushing et al.ll2005l ; iRayner et al.ll2009l ). We find that the candidate is a M4.5 dwarf based 
on this spectrum. Thus, coupled with the measured magnitudes and colors (see Table |2j), it 
is consistent with being a bona fide companion. 
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Fig. 9. — SpeX prism spectrum of the candidate companion to HIP 97944 (black solid 
l ine) compared to M3-M6 dwarf spectra (red dotted lines) from the IRTF Spectral Library 
(ICushing et al.ll2005l ; iRayner et al.ll2009l ). We find that the candidate is a M5 dwarf based 
on this spectrum. Thus, coupled with the measured magnitudes and colors (see Table |2j), it 
is likely a bona fide companion. 



31 




Fig. 10. — J-band magnitude limits as a function of radial separation from the central 
spectroscopic binaries for the fields of HIP 39064 (left) and HIP 44248 (right), with 50% 
completeness displayed as the solid red line and 90% completeness as the dashed green line. 
HIP 39064 (V = 7.7 mag near the median value of primary magnitudes for our survey), 
represents our typical sensitivity, whereas HIP 44248 is at the bright end of our sample, with 
V = 4.0 mag. Note that, as expected, the sensitivity as a function of separation is degraded 
for the brighter star. 
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Fig. 11. — Histogram of 50% and 90% completeness J band limiting magnitudes for all 
observed fields. There is a wide range of completeness limits, due to the changes in seeing 
and primary brightness across all the observed fields. The median J band completeness limit 
at 50% is 19.1 mag and at 90% is 18.2 mag. 
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Fig. 12. — Histogram of the period of the central spectroscopic binaries (in Log days) versus 
the calculated companion frequency. Note the large spike in com panion frequency at the 
smal lest binary periods. This is c onsistent with the predictions of (jSterzik fc Durisenl 120031 ) 
and (jDelgado-Donate et al.ll2004J ). 
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Fig. 13. — Scatter plot of the period of the central spectroscopic binaries with tertiaries 
against the projected separation of its tertiary companion. Note that there is little visual 
evidence of a correlation between them. 
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K5 


03 


May 


2007 


KPNO 


77725 


15 


52 


08 


.24 


+10 52 28.1 


44. 


.27 


± 


1 


.56 


-267.30 ±2.17 


-238.36 ± 1.76 


1014 


,500 


9 


.38 


MOV 


06 


June 


2007 


CTIO 


77801 


15 


53 


12. 


.19 


+13 11 52.8 


57. 


.27 


± 





.88 


-150.56 ± 1.32 


-562.69 ±0.98 


138. 


.603 


6 


.10 


G0IV 


07 


June 


2007 


CTIO 


78527 


16 


01 


53 


.70 


+58 33 52.0 


47. 


.79 


± 





.54 


-320.07 ±0.50 


334.96 ±0.53 


3. 


071 


4 


.01 


F8IV-V 


06 


May 


2007 


KPNO 


78709 


16 


04 


03. 


.71 


+25 15 17.4 


46. 


.56 


± 





.89 


-488.79 ±0.58 


696.64 ±0.73 


4450. 


.800 


7. 


.06 


G8V 


04 


May 


2007 


KPNO 


79607 


16 


14 


41 


.04 


+33 51 31.8 


46. 


11 


± 





.98 


-266.47 ±0.86 


-86.88 ± 1.12 


1. 


.140 


5 


.64 


GOV 


06 


May 


2007 


KPNO 


80686 


16 


28 


27. 


.80 


-70 05 04.8 


82. 


.61 


± 





.57 


199.89 ±0.31 


110.77 ±0.51 


12. 


.976 


4 


.90 


GOV 


05 


June 


2007 


CTIO 


80925 


16 


31 


30. 


.35 


-39 00 41.3 


40. 


.60 


± 


1 


.75 


-428.05 ± 1.47 


-333.41 ± 1.43 


31. 


.846 


7. 


.25 


dKl 


04 


June 


2007 


CTIO 


81693 


16 


41 


17 


.48 


+31 36 06.8 


92. 


.63 


± 





.60 


-462.58 ±0.59 


345.05 ±0.66 


12596 


.100 


2. 


.80 


G0IV 


04 


May 


2007 


KPNO 


82860 


16 


56 


01 


36 


+65 08 04.8 


66 


.28 


± 





.48 


238.05 ±0.38 


50.84 ±0.57 


52. 


.109 


4 


.88 


F6V 


03 


May 


2007 


KPNO 


83608 


17 


05 


20. 


.12 


+54 28 12.2 


37. 


.08 


± 





.89 


-66.00 ±0.98 


73.86 ± 1.03 


42309 


.000 


5 


.83 


F7V 


04 


May 


2007 


KPNO 


85667 


17 


30 


23. 


.87 


-01 03 45.0 


60. 


.80 


± 


1 


.42 


-126.64 ± 1.72 


-172.00 ±0.91 


16830 


.400 


5 


.30 


G8IV-V 


04 


June 


2007 


CTIO 


86036 


17 


34 


59 


.59 


+61 52 28.4 


70. 


.98 


± 





.55 


277.38 ±0.54 


-525.62 ±0.60 


27087 


.000 


5 


.23 


G1V 


06 


May 


2007 


KPNO 
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Hipparcos RA 


DEC 


7T a 




M<5 a 


Pcriod b V b SpT c UT Date Obs. Telescope 


Number (h am as) 


(d m s) 


(mas) 


(mas/yr) 


(mas/yr) 


(days) (mag) 



86901 


1 7 


36 

OU 


r i7 


OQ 

■ U -J 


+68 45 95 9 


49 69 + 53 

4zZ(.Uz-i —1— U.OO 


1 Q4 _|_ n r.q 


Q21 05 + 69 


K 

O ■ 


980 


4 


80 

■ OO 


F5V 

X O V 










86400 

OU4UU 


1 7 


3Q 

O -J 


1 fi 


91 


+03 33 1 8 9 

\^\JO OO J-O.J7 


93 36 + 1 25 


— 1 79 67 + 75 

J- 1 iJ .U I —1— O. 1 O 


_qc 24 + 59 


83 

OO. 


798 


fi 
\j 


59 

. Oi 


K3V 

-LVO V 


06 

OO 


Tnnn 

O UIll_< 


9007 

ZjUO 1 


CTTO 


8fi799 

OU ( 


1 7 


4H 


1 (J 


fi4 


-1-91 ^fi ^ 9 


AO A ^ 4- n QR 


J_zjO.-LO ZIZ l.UU 


— fil Q 84 + RR 


9^8 


400 


7 


^0 

. OU 


KOV 

IV U V 


07 
u / 


J line 


9007 

ZiUU 1 




OOUU -L 


1 8 


u o 


97 


.21 


+09 Qf) OS 8 


1 96 69 + 1 38 

_L£7U.UZi I -L.OO 


1 94 Kfi _i_ 1 it; 


—969 66 + 91 

iyU^j.OO —1— 0.t7-L 


^91 


R01 


4 


09 


KOV 

.LVO V 


Oft 
00 


Tl lTlO 
O LIIHj 


9007 

ZjOO 1 


CTTO 


89937 


1 8 


21 


V J £i 


34 


+79 44 01 3 


104 11 + f) 48 


wi 08 + 49 

OO-L.OO —1— \J .^±cr 


_qci cq _i_ 4fi 

OO-L.OC/ —1— O.TrO 


980 


550 

OOO 


O 


57 

■ O I 


F7V 

J. I V 










90355 

£70000 


18 


26 


10 


.00 


+08 46 39 

I OO i U Ot/.O 


O I . U t: _l_ IiIU 


104 «q + 1 45 


_4ko 10+1 97 

t:(JU. lu —1— J. . — 1 


293 


500 


7 


.78 


G7V 


05 


June 


2007 


CTIO 


91009 


18 


33 


55 


.60 


+51 43 1 1 7 

I O -L i O J. J. . 1 


60 90 + 73 

UU.i/O —1— O . 1 o 


iqc «o + 71 


—394 89 + 77 

OZi4:.Oi7 —1— O. 1 1 


5. 


.975 


g 


.04 


dMOe 










92919 


18 


55 


53 


.14 


+23 33 26 4 


4« «4 -f i hq 

t:U.ui —1— J-.VJO 


1 30 79 + 59 

-LOO. 1 U —1— O . O U 


—983 07 + 74 


2, 


.879 


8. 


,30 


KOV 


07 


June 


2007 


CTIO 
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18 


57 


01 


.61 


+32 54 04 6 

I OZj Ot: Ot: . O 


rr 7« 4. n 54 


2Q2 ok + 46 


_i4Q q7 _|_ 4q 


22423 


.000 


5. 


.22 


GOV 










931 74 

i?0 J. 1 4: 


1 8 
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(JO 


40 


47 
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O / OU ZiO . o 
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00.4:0 -I— \J .ij Li 


_i Q2 95 + 1 

lOZi.^U —1— X.OO 


_i 1 45 + 77 




O ■ 


591 

Oi7 X 


4 


74 


F0V 

1 U V 


04 

04 


T 

j unc 


9007 

iOU i 


CTTO 

\J _L 1 VJ 


93996 




07 

U / 


32 

Oi 


39 


+30 15 16 2 

n^oo j.o 


35 70 + 78 
oo. i u -i— o. i o 


111 or + 63 

III.l/U _i_ 0.00 


1 0*^ 03 + 68 

-LOO.OO —1— O.OO 


2. 


1 31 


7 


63 

-UO 


O8V 

OO V 








00 




I Q 


U ( 


O i 


9R 


4-1 ^.1 14 Q 
I ± u o± i^.y 


47 79 _i_ n 77 




_^04 4fi 4- fiD 
OU'i.'lU zc u.uu 


91 
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a 
u 


OS 

.UO 


vi -J V 
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uu 
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19 


26 


25 


.98 


+49 27 55 1 


qq 70 4. 1 no 

Ot/. 1 O —1— l.UO 


457 1 q + 1 09 

t:0 1 . J. t/ —1— l.UJ 


71 Q Q7 _L 1 14 


166. 


.360 


8 


.01 


K3V 










95995 






(JO 


54 


+ 58 35 13 1 

T<JO OO J.O.-L 


Oi/.O^t _l— u.u^ 


— 51 04 + fiS 


— ^07 ^4 4- fiR 

Oc7 1 .04: _l_ O.OO 


4Q1 


900 


u 


59 


K1V 










070/1 /I 

y fy44 


1 O 


04 


1 1 


.OZ 


OQ C/? O/l Q 

-ZO 00 z4.o 


U.o4 zt U.ol 


— Izz.D ( i U. 1o 


/inn j- n e:q 
— 4Uy.oD zt U.oo 


40. 


,ol ( 


D 


ID 


1VO V 


n 1 
U4 


June 


onnv 
zUU ( 




98416 


19 


59 


47 


.49 


-09 57 26.2 


40.75 + 1.35 


— 246.73 + 2.31 


-392.36+ 1.63 


1786 


.270 


6 


.22 


F8V 


06 


June 


2007 


CTIO 


101382 


20 


32 


51 


.76 


+41 53 50.6 


44.99 ± 0.64 


156.89 + 0.53 


452.80 + 0.47 


57. 


.321 


7. 


.09 


G9V 










101750 


20 


37 


20 


.82 


+75 35 56.7 


36.16 + 0.97 


309.20+ 1.12 


539.49 + 0.97 


0. 


278 


7. 


.23 


KOV 










101769 


20 


37 


32 


.87 


+14 35 42.7 


33.49 ± 0.88 


118.28+1.33 


-47.65 + 1.09 


9733 


.700 


3 


.62 


F5IV 


07 


June 


2007 


CTIO 


101955 


20 


39 


37 


.20 


+04 58 18.7 


53.82 + 2.21 


862.35 + 2.15 


67.57+ 1.78 


920. 


.200 


7. 


.89 


K5V 


05 


June 


2007 


CTIO 


102431 


20 


45 


21 


.12 


+57 34 47.0 


36.87 + 0.46 


-62.95 + 0.41 


-235.56 + 0.45 


523 


.360 


4 


.51 


F8IV 










103655 


21 


00 


05 


.35 


+40 04 13.0 


66.21 + 2.54 


614.37 + 2.07 


-247.16 + 2.45 


10777 


.100 


10 


.10 


dM3e 










104858 


21 


14 


28 


.79 


+10 00 27.8 


54.11 + 0.85 


42.32 + 0.89 


-303.43 + 0.57 


2082. 


.100 


4 


.49 


F5V 


05 


June 


2007 


CTIO 


105312 


21 


19 


46 


.00 


-26 21 07.2 


53.40 + 1.09 


-582.35 + 1.11 


-357.67 + 0.62 


21. 


.346 


6 


.55 


G5V 


04 


June 


2007 


CTIO 


105406 


21 


21 


01 


.44 


+40 20 44.1 


37.64 + 0.59 


-18.89 + 0.43 


-208.92 + 0.47 


3. 


243 


6 


.43 


F8V 










107089 


21 


41 


28 


.47 


-77 23 22.1 


47.22 + 2.90 


64.83 + 2.63 


-240.38 + 2.12 


1020. 


.000 


3 


.75 


K0III 


06 


June 


2007 


CTIO 


107556 


21 


47 


02 


.29 


-16 07 35.6 


84.58 + 0.88 


263.26+ 1.23 


-296.23 + 0.67 


1. 


023 


2 


.83 


Am 


04 


June 


2007 


CTIO 
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Hipparcos 




RA 


DEC 


7T a 






Period 13 


V b 


SpT c 


UT Date Obs. 


Telescope 


Number 




am as) 


(d m s) 


(mas) 


(mas/yr) 


(mas/yr) 


(days) 


(mag) 








109176 


22 


07 00.47 


+25 20 42.2 


85.06 + 0.71 


296.73 + 0.56 


26.93 + 0.69 


10.213 


3.76 


F5V 






111170 


22 


31 18.22 


-06 33 17.6 


39.18 ± 1.83 


160.65 ± 1.15 


-108.63 ± 0.80 


630.140 


6.27 


F7V 






111802 


22 


38 45.29 


-20 37 15.4 


115.71 + 1.50 


450.58+ 1.61 


-79.86+ 1.15 


4.083 


9.09 


dMle 


07 June 2007 


CTIO 


113718 


23 


01 51.54 


-03 50 55.4 


59.04 + 3.40 


395.50+ 1.21 


-207.11 + 1.17 


454.660 


7.46 


K4V 






113860 


23 


03 29.76 


-34 44 58.6 


34.98 + 0.80 


74.80 + 0.81 


84.45 + 0.70 


178.318 


5.10 


F0V 


07 June 2007 


CTIO 


114379 


23 


09 57.23 


+47 57 30.0 


39.56 + 7.67 


147.06 + 5.75 


12.42 + 6.72 


3.033 


7.91 


K2V 






115126 


23 


19 06.51 


-13 27 30.4 


48.22 + 5.25 


265.64 + 6.17 


-44.87+5.22 


2323.600 


5.21 


G5IV 






116584 


23 


37 33.71 


+46 27 33.0 


38.74 + 0.68 


159.22 + 0.33 


-421.46 + 0.51 


20.521 


3.88 


G8III-IV 






116727 


23 


39 20.85 


+77 37 56.2 


72.50 + 0.52 


-48.85 + 0.48 


127.18 + 0.44 


24135.000 


3.23 


K1IV 






117712 


23 


52 24.52 


+75 32 40.2 


92.68 + 0.55 


341.82 + 0.53 


41.88 + 0.47 


7.753 


6.40 


K3V 







CO 
CO 



a Data in these columns derived from the Hipparcos database (jPerrvman et al.lll997l ). 

b Data in these columns derived from the Ninth Catalog of Spectroscopic Binaries (|Pourbaix et al.ll2004l) . 



= Spectral Types taken from those given in the Ninth Catalog of Spectroscopic Binaries (jPourbaix et al.l 12004 ) 



Table 2. Candidate and Known CPM Companions 



Primary 


2MASS RA 


2MASS DEC 






Sep. 






J - K s c 




J - H 


C 


H - K s c 


SpT d 


SpT 


Notes 


Comps 


Hip. # 


Terti £Lry 




Tertiary 


(mas /yr) 


(mas/yr) 


(AU) 


(mag) 




(mag) 




(mag) 




(mag) 


Tert. 


Pri. 




Stati 


13081 


02h48m09 


.7s 


+27d04m25s 


237.1 ±40 


-122.1 ±44 


160 


8.980 ±0.022 





.862 ± 0.031 


0. 


.554 ±0. 


.040 


0.308 ± 0.040 


M4.5 


K1V 


l c 


Ye 


16846 


03h36m46 


,8s 


+00d35ml5s 


-28.2 ± 32 


-131.6 ±32 


220 
















K6 


K1V 


5 f 


Ye: 


21482 


04h36m44 


.9s 


+27d09m51s 


254.1 ±44 


-201.7 ±40 


2180 


13.340 ± 0.038 





.462 ± 0.079 


0. 


.366 ±0. 


069 


0.096 ± 0.090 


DC 


dK5pe 


is 


Ye: 


36850 


07h34m37 


.4s 


+31d52ml0s 


-207.6 ±28 


-96.0 ± 28 


1130 


5.079 ±0.018 





.837 ±0.027 





.653 ±0. 


028 


0.184 ±0.029 


M0.5 


A1V 


5 h 


Ye: 


41211 


08h24m33 


,8s 


-03d44m34s 


-211.5 ±36 


-14.2 ± 36 


980 


9.267 ±0.027 





.841 ± 0.037 





.547 ±0. 


035 


0.294 ±0.034 


M4.5 


F1V 


2 f 


Ye: 


41211 


08h24m52 


.3s 


-03d41m02s 


-217.0 ±36 


-17.0 ± 36 


9640 


9.374 ±0.026 





.918 ±0.032 





.566 ±0. 


033 


0.352 ±0.028 


M5.5 


F1V 


1' 


Ye: 


42172 


08h35m51 


3s 


+06d37m22s 


-120.7 ±48 


-133.3 ±44 


230 


4.033 ± 0.024 





.380 ±0.032 


0. 


.170 ±0. 


.036 


0.110 ±0.034 


G5 


F5V 


5J 


Ye: 


45170 


09hl2ml4 


.7s 


+14d59m40s 


-565.5 ±32 


279.3 ± 40 


860 


13.955 ±0.078 


1 


.469 ±0.100 





.891 ±0. 


122 


0.578 ±0.103 


L8 


G9V 


l k 


Ye: 


46509 


09h29m09 


.2s 


-02d45m03s 


138.7 ±36 


-20.6 ± 40 


1120 


4.529 ±0.029 





.557 ± 0.036 





.391 ±0. 


.057 


0.166 ±0.054 


K0 


F6V 


5J 


Ye: 


72603 


14h50m27 


.4s 


-16d03m23s 


-112.5 ±28 


-105.6 ±32 




14.603 ±0.120 


1 


.065 ± 0.257 


0. 


.535 ±0. 


206 


0.530 ±0.281 




F3V 


3 


Nc 


75312 


15h23m22 


.6s 


+30dl4m56s 


133.1 ±32 


-183.2 ±28 


3640 


14.706 ± 0.099 


1 


.708 ±0.120 


1 


.128 ±0. 


.128 


0.580 ±0.105 


L8 


GOV 


I 1 


Ye: 


75718 


15h28ml2 


.2s 


-09d21m28s 


83.5 ± 28 


-355.9 ± 28 


980 


4.504 ±0.027 





.533 ± 0.034 





.443 ± 0. 


.048 


0.090 ± 0.045 


K2 


K1V 


5S 


Ye: 


83608 


17h05m20 


3s 


+54d28m00s 


-74.5 ± 24 


89.5 ±36 


340 


6.866 ± 0.043 





.591 ± 0.046 





.356 ±0. 


.051 


0.235 ±0.032 


?? 


F7V 


2 m 


Ye: 


86722 


17h43ml5 


.4s 


+21d36mlls 


-98.7 ± 40 


-631.5 ±40 


530 


9.650 ±0.025 





.811 ± 0.031 





.495 ± 0. 


033 


0.316 ±0.028 


M4.5 


KOV 


3 


Ye: 


92919 


18h55m50. 


,4s 


+23d36m51s 


-29.0 ± 28 


-154.9 ±28 




10.797 ±0.021 





.816 ± 0.029 





.579 ±0. 


.030 


0.237 ±0.029 




Kofe 


4 


Nc 


97944 


19h54m20 


6s 


-23d56m40s 


-151.7 ±32 


-400.8 ± 32 


590 


9.027 ±0.021 





.952 ± 0.031 





.598 ±0. 


032 


0.354 ± 0.033 


M5 


K3Y 


3 


Yei 


101769 


20h37m33 


,7s 


+14d35m32s 


50.9 ±36 


-55.6 ±36 




8.086 ± 0.046 


-0.008 ±0.096 





.358 ±0. 


.070 


-0.366 ± 0.099 


GO 


F5IV 


3 


Nc 


111802 


22h38m45 


.3s 


-20d36m52s 


443.1 ±40 


-25.2 ±36 


200 


7.661 ±0.024 


0.853 ±0.029 





.527 ±0. 


.062 


0.326 ±0.059 


M3.5 


M1.5V 


l e 


Ye: 



a isted proper motion measurements are derived from the data presented in this paper save for those objects with a [5] in the follow-up column. Those values were taken from Hipparcos, as 
objects were saturated in our data. 

b Absolute J magnitudes derived from 2MASS photometry and Hipparcos parallaxes. 

C A11 photometry listed is from the 2MASS PSC. 

d Companion spectral types were based either on SpeX prism spectroscopy obtained for this work (objects with [2] or [3] in the Notes column), or were already known in which case the spectrE 
listed is from the paper referenced by the lettered superscripts in the Notes column (objects with [1] or [5] in the Notes column). 

jReid et al.l Jl995l) 

iTokovinin et al.l J2006h 

^Makarov et al.1 feOOSl) 

^Joy fc Sanfordl jl92d) 

jReid et al.l J2003I) 



] Jl99Gh 
(2001) 

:et fc Nyj j2002h 

) Previously known companion, no additional follow-up data required. (2) Previously known candidate companion, no Spectral Type, SpeX data acquired. (3) No previous data, SpeX 
. (4) Additional photometry acquired from Vizier database, not a companion. (5) Previously known bright companion. Saturated in these data, found via literature search using the 
ibasc. 
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Table 3. Companion Fraction by Spectral Type 



Spectral 
Type 


Number of 
Primaries 


Number with 
Tertiaries 


Companion 
Fraction 


B 

A 
F 
G 
K 

M 


1 

4 
23 
21 
24 
4 



1 
5 
2 
6 
1 


0.0 + gV% 
25.0+?^% 
21.7+1°?% 
9.511°?% 
25.0+l° 6 5 % 
25.0±?S;|% 



Table 4. Companion Fraction by Spectroscopic Binary Period 



Period Number of Number with Companion 
Log Days Primaries Tertiaries Fraction 



-0.5 


2 





0.0tf 8 % 


0.5 


12 


6 


50.0 + ^|% 


1.5 


19 


1 


5.3tf 7 2 % 


2.5 


18 


2 


n.i+iy^o 


3.5 


19 


6 


31.6™% 


4.5 


7 


2 


28.6+f 3 7 % 



